The stability of strained H:Si(105) and H:Ge(105) surfaces 
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We report atomic scale studies of the effect of applied strain and hydrogen environment on the 
reconstructions of the (105) Si and Ge surfaces. Surface energy calculations for monohydride- 
terminated (001) and (105) reconstructions reveal that the recently established single- height re- 
bonded model is unstable not only with respect to (001), but also in comparison to other mono- 
hydride (105) structures. This finding persists for both Si and Ge, for applied biaxial strains from 
-4% to 4%, and for nearly the entire relevant domain of the chemical potential of hydrogen, thus 
providing an explanation for the recently observed H-induced destabilization of the Ge(105) surface. 



The epitaxial system Ge/Si(001) has been the focus 
of intense investigations for more than two decades, act- 
ing both as a test bed for our fundamental understand- 
ing of the strained-layer growth, as well as a technolog- 
ical launching pad for promising optoelectronic devices 
based on Ge/Si quantum-dots. While the understand- 
ing of the formation of quantum dots (islands) has pro- 
gressed rapidly 1], the desire to further diversify and con- 
trol growth morphologies has triggered studies of Ge/Si 
epitaxy in the presence of other species. In particular, 
hydrogen was shown to have a surfactant effect on the 
deposition of Ge on low-index silicon surfaces: the island 
formation was suppressed in the presence of atomic hy- 
drogen, with the growth switching to the layer-by-layer 
mode 0,0- In the H-mediated Ge/Si(001) epitaxy, Si 
atoms tend to segregate at the surface, and their ex- 
change with Ge atoms is reversible upon hydrogen des- 
orption Q . Recently, experiments have broken ground in 
a different direction to address the influence of hydrogen 
on high- index epitaxial systems. Fujikawa and cowork- 
ers used scanning tunneling microscopy and electron en- 
ergy loss spectroscopy to investigate the H adsorption on 
Ge/Si(105), and demonstrated the destabilizing effect of 
hydrogen on the surface p|. 

Motivated by these compelling experiments 0, we 
have studied the influence of applied biaxial strain and 
chemical potential of hydrogen on the surface energy of 
(105) and (001) reconstructions of Si and Ge. Our model 
potential calculations predict that for a wide range of ap- 
plied strain and H chemical potential, the single-height 
rebonded |E IS El structure with monohydride termi- 
nation has higher surface energy than the H-terminated 
(001). Interestingly, this rebonded-step structure (RS 
or SR 0) also becomes unstable with respect to the 
single-height unrebonded (SU) model that was originally 
proposed for the configuration of (105) hut facets [l0(- 

In addition to SU and SR, we have studied all other 
(105) structures although, for brevity, we report 

results only for the SU, SR, DT, and DU models (de- 
picted in Fig.nj. Using the Tersoff potential param- 
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FIG. 1: (Color online) Monohydride (105) reconstructions 
with two different densities of hydrogenated bonds (Hb), 
1.57 Kb I a 2 (SR, DT) and 2.35Hb/a 2 (SU, DU), where a is 
the bulk lattice constant of Si or Ge. 



eterized for the Si-H and Ge-H systems [3, [^| , we have 
performed full relaxations of 150 A-thick periodic slabs 
with either clean or H-passivated reconstructions under 
applied biaxial strains e in the range —4% < e < 4%. 
With the dependencies on strain (e) and H chemical po- 
tential made explicit, the surface energy 7 of a hy- 
drogenated reconstruction has been calculated as 

7 (e, li h ) = (E(e) - Nfi(e) - N H » H )/A(e), (1) 

where E is the total energy of the slab of area A, N 
is the number of Si or Ge atoms, Njj is the number of 
H atoms passivating the surface, and \i is the chemical 
potential (bulk cohesion energy) of the Si or Ge atoms. 
Guided by experiments 5] , we have only considered here 
monohydride terminations and estimated the range \±h 
for which the surface hydrogenation becomes thermody- 
namically favorable. We found that the lowest value hh 
at which an H-passivated surface becomes favorable over 
the pristine surface is only weakly dependent on the re- 
construction and strain, so we use one chemical potential 
range for each material: —3.1 eV < hh < 0.0 eV for Si 
and -2.8 eV < fi H < 0.0 eV for Ge Efl. 

Our main findings are summarized in Fig. which 
shows the energies of Si and Ge surfaces as functions 
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FIG. 2: (Color online) Surface energies (meV/A 2 ) of clean 
(a,b) and monohydride (c,d) Si and Ge surfaces plotted ver- 
sus biaxial strain e (%). The bottom panels show surface en- 
ergy curves at the upper and lower bounds of the H chemical 
potential range for Si (c) and for Ge (d). 



of the biaxial strain e. The results for clean sur- 



faces [Fig . Efa,b)] are consistent with previous studies 
El El El which point out that under compressive strain 
the SR reconstruction becomes more stable than the fiat 
(001) surface. In the case of hydrogenated Si and Ge 
surfaces we find that: 

(a) the ordering of surface energies corresponding to 
SR-(105) and (001)-2 x 1 is reversed upon hydrogenation 
for the entire range of pu and e, for both Si and Ge. 

(b) different (105) reconstructions have surface ener- 
gies whose plots vs. e are grouped (bunched) according 
to the number of hydrogenated bonds per area (e.g., SR 
and DT, SU and DU). 



TABLE I: Surface energy ratios of monohydride reconstruc- 
tions at e = and /jh=0.0 eV, compared to the ratios of 
densities p of passivated bonds. 

SU/(001) (001)/SR SU/SR SR/DT 
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1.24 


1.65 


2.04 


1.07 
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ratio, Ge 


1.26 


1.75 


2.20 


1.09 


P 


ratio 


1.18 


1.27 


1.50 


1.00 



(c) the slope dj/de is positive for all monohydride re- 
constructions, irrespective of what its sign was for clean 
surfaces. 

On general grounds, observation (a) is consistent with 
ab initio calculations which estimate that the average 
binding energy of H on Ge/Si(105) is 0.1 eV smaller than 
that of H on Ge(001) 5]. Still, this estimate may not 
translate readily into a surface energy difference (or ra- 
tio) between the monohydride terminated SR and (001). 
We therefore pursue further the origins of the reversal of 
energetic ordering of SR and (001) upon hydrogenation. 
The ab-initio binding energy of the Ge-H bond is 2.83 eV 
[TEj . somewhat smaller than that of the Si-H bond. Since 
this binding energy is an order of magnitude larger than 
the difference Q in average binding of H on Ge/Si(105) 
and H on Ge(001), the key factor in determining the en- 
ergetic ordering of monohydride reconstructions is the 
number of passivated bonds per area. Indeed, if we focus 
on surface energies at e = 0% and hh = 0.0 eV, we find 
that their relative ratios can be largely accounted for by 
the ratios of the areal density p = Nh/A of monohydride 
bonds, as shown in Table |I] for both Si and Ge. Thus, 
the SR structure becomes unstable with respect to (001) 
mainly because it has a lower density of H-bonds (Hb) at 
the surface (i.e., 1.57 Hb/ a 2 as opposed to 2.00Hb/a 2 for 
(001), where a is the bulk lattice constant of either Si or 
Ge). 

Continuing the discussion at pn — 0.0 eV, when the 
H-bond density p is the same for two different reconstruc- 
tions, then their energetic ordering is dictated by subtler 
differences in surface stress and in adsorption (binding) 
energy of H. In general, these factors have less influence 
on surface energies than the H-bond density (Table 
which is why, e.g., SR and DT curves are bunched to- 
gether and are relatively far from the other surface en- 
ergy curves at pn = 0.0 eV (Fig. EJ). We have found 
this bunching trend [observation (b) above] for all the 
(105) reconstructions that have the same density of 
passivated bonds. 

Variations of pn in the stated ranges lead, via Eq. (QJ, 
to variations of the energetic separation between same-p 
bunches, but leave the surface energy gaps within each 
bunch unchanged (e.g., the gap between SR and DT does 
not depend on pn)- Because the lower bounds of pn 
(Fig. [2J) necessarily have magnitudes that are similar to 
the Si-H or Ge-H binding energies, the separation be- 
tween equal- p bunches can be wiped out in those limits 
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[refer, e.g., to the curves corresponding to hh = —3.1 
eV in Fig. Efc)]. On the other hand, an increase of \±n 
by 0.3-0.5 eV above the lower bounds of fin is sufficient 
for the bunches to be clearly distinguishable, thus for the 
density of passivated bonds to set in as the dominant 
factor for surface energetics. 

We now turn to discussing the observation (c) listed 
above. The positive slope of the surface energies in 
Fig. Etc, d) is directly related to an increase in the diago- 
nal components of the surface stress of upon hydrogena- 
tion. We have calculated the surface stress components 
for all clean and monohydride Si (001) and (105) models 
and found that each of the two diagonal components an 
and (J22 increases by as much as 40.0-65.0meV/A 2 upon 
hydrogenation. In the case of Ge surfaces, the trend is 
the same but the increases are smaller, in the range of 
15.0-25.0 meV/A 2 . Noting that stress components do 
not depend on the increase of an and 022 is caused 
by the stretching of surface bonds. Indeed, we have ver- 
ified for Si surfaces that the length of the dimer bonds 
is stretched by about 2.5% upon passivation, while the 
bridge-bonds made by the rebonded atoms is increased 
by about 1.6%. On Ge surfaces, dimers are stretched by 
2.0% and bridges by 0.8% , which explains the smaller in- 
crease of stress components calculated for Ge. Ab initio 
calculations on H:Ge/Si(105) indicate that the bridge- 
bonds actually decrease by 1.5% 5] for the Ge/Si(105) 
system, which is also consistent with the smaller stress 
variations for Ge surfaces compared to Si ones. 

Before concluding, we comment as to why the empiri- 
cal potentials can capture the main energetic trends for 
hydrogenated surfaces. We note that Tersoff pjj po- 
tentials have already been found to predict the correct 
lowest-energy reconstruction and strain behavior in the 
absence of hydrogen 0, 0, [HI , although most surface 
energy gaps were overestimated. For passivated surfaces, 
the hybridization state of all atoms is sp 3 5], and thus 
major departures from the fitting databases of these po- 



tentials do not occur. Subtle electronic effects are there- 
fore unlikely to dominate the energy ordering of vari- 
ous monohydride reconstructions, which renders them 
tractable at the level of Tersoff potentials 0, 0, Il5| . 
We hope that the results obtained here will be followed 
by further studies at the ab initio level, and we expect 
quantitative improvements to emerge mainly for compar- 
isons between clean and H-passivated surfaces. 

In conclusion, we investigated the relative stability of 
strained (105) reconstructions with monohydride termi- 
nations. On purely thermodynamic grounds, our calcu- 
lations suggest that the monohydride H:Ge/Si(105)-SR 
structure (5| is unstable for any Ge or Si-Ge coverage (or 
applied strain), not only for coverages below a threshold 
value. However, the kinetic barriers against SR decaying 
into other monohydride reconstructions and the compres- 
sion in the epitaxial layer make the H:Ge/Si(105) experi- 
mentally observable as reported in Ref. p . Since the SU 
model has the largest possible surface density of mono- 
hydrides and a structure based on single- height steps, 
it has been found here to also have the lowest surface 
energy for the relevant ranges of strain and H chemical 
potential. Thus, it is conceivable that the SU model may 
be confirmed experimentally by pursuing different surface 
preparation recipes that diminish the influence of kinetics 
while avoiding the formation of dihydrides. The results 
presented here could also be relevant for future experi- 
ments on the Si(105) surface, which has long been known 
to be rough and disordered on atomic scales 0, IT^ . 
If maximal monohydride coverage can be achieved on 
Si(105), the present results indicate that Si(105) would 
become ordered in the presence of H, with a reconstruc- 
tion given by the SU model. 
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